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Abstract. Despite the rapid observational progress in the study of young massive star 
clusters, the formation of globular clusters still remains poorly understood. Yet, it is 
emerging that globular cluster formation is intimately linked to the formation of the 
Galaxy. I discuss a generic scenario of the cluster formation within progenitor galaxies, 
based on the available observational constraints. The oldest clusters formed around 
redshift z = 7, but the process continued at least until z = 3. Because of their high 
density, globular clusters survived when their progenitor galaxies were disrupted by 
the Galactic tidal field. 

1 Observational Constraints 

Building a successful physical model of globular cluster formation has proved to 
be an elusive topic, as revealed by the multitude of effort and the lack of result 
in the last 30 years. In essence, this is the point where the large-scale galaxy 
formation problem meets the small-scale star formation problem. In order to 
try to overcome the uncertainties associated with both, we should seek as many 
clues as possible from observations. 

Fortunately, recent observational studies of star-forming regions in the Galaxy 
and young star clusters in interacting galaxies have improved our understand- 
ing of star formation. The majority of stars forms in clusters and associations 
of various sizes. The hierarchy ranges from the young OB associations to the 
old massive globular clusters, with no special scale between 10 and 1O 6 M : 
dN/dM oc M~ a , a = 1.5 — 2 [pHl9|]. Thus globular clusters seem to represent 
the most massive end of the general hierarchy. 

The observed luminosity function (LF) differs from the initial one due to 
the dynamical evolution. Small- mass clusters (M < 10 5 M Q ) are gradually de- 
stroyed by stellar two-body relaxation and tidal interactions ]l8|,p|,p^H^] , while 
very massive clusters sink to the center via dynamical friction. However, most of 
the high-mass clusters are essentially unaffected by the evolution, and therefore, 
preserve the shape of the initial LF. 

Observations and modeling of old Galactic and extra-galactic globular clus- 
ters reveal similar general properties that depend little on the vastly different 
environment of the host galaxies (small and large, spiral and elliptical): 

• Age. The oldest globular clusters are old, but realistically old: tec ~ 12.5 
Gyr, with the 95% probability range 10.2-15.9 Gyr fuj. There is no problem 
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now with the stars being older than the Universe, and their age can be used 
to place constraints on the epoch of globular clusters formation: 

tf = tff — tec, 

where tu ~ 13.5 Gyr is the Hubble time for the concordance cosmological 
model [pj. With all the associated uncertainties, oldest clusters must have 
formed between redshifts Zf = 7 and 3 

• Metallicity. A large spread of metallicity among old clusters, from 0.01 to 
0.3 Zq, with a median ~ 0.03 Zq. The material from which globular clusters 
formed must have been pre-enriched by heavy elements. A narrow range of 
metallicity within individual clusters (<5[Fe/H] < 0.1) indicates that each 
globular cluster formed in a single burst. 

• Distribution. Spherically symmetric distribution and little rotation of the 
Galactic globular cluster system point to a different origin than most stars 
in the rotationally-supported disk. Also, the efficiency of globular cluster 
formation is universally low [|l3) : the ratio of the mass of globular clusters to 
the total baryonic (stellar + gaseous) mass of their host galaxy is the same 
in giant and dwarf ellipticals, €qc = Mcc/M^ ss 0.0026 ± 0.0005. 

2 A generic hierarchical scenario 

Early models of globular cluster formation |lj],|]|| focused on local physical pro- 
cesses, assuming a static environment of the present Galaxy. In the last decade 
the paradigm of galaxy formation has changed significantly. We now have good 
evidence that galactic halos form as a result of gravitational growth and inter- 
action of primordial fluctuations within the cold dark matter model . Small 
objects collapse first and merge into larger systems, extending the hierarchy to 
progressively higher masses. The present halo of the Milky Way formed in dozens 
of mergers of smaller progenitors . 

At redshift Zf ~ 7, globular clusters must form within small progenitor galax- 
ies, because large galaxies do not exist yet. The primordial gas would cool by the 
hydrogen recombination line emission within the first halos with the virial tem- 
perature T v ; r > 10 4 K and settle into a rotationally-supported, self-gravitating 
disk. Local instabilities would break the disk into dense clumps, which could 
further fragment into stars. Star formation in the local universe proceeds in thin 
disks, and it is natural to assume the same for the high-redshift progenitors. 
The overall spherical geometry of the globular cluster system follows from the 
subsequent mergers and tidal disruption of the progenitor galaxies. 

For illustration, consider a simple scenario of the first episode of star forma- 
tion in a progenitor galaxy. Assume that the scale-free mass function of giant 
molecular clouds translates into a similar cluster mass function, dN c i/dM oc 
M~ a . For a < 2, the total mass of the star cluster hierarchy within a progenitor is 
dominated by the most massive cluster, M max . Combined with the observed for- 
mation efficiency cgc, the total baryon mass of the progenitor could be estimated 
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Fig. 1. Cumulative number of progenitors of the present Milky Way halo more massive 
than M at three different redshifts. Overplotted by asterisks is the luminosity function 
of the massive Galactic globular clusters rescaled by a factor (fb^Gc)' 1 ~ 3000, for 
a constant mass-to-light ratio M/Ly = 2. If some progenitors had higher formation 
efficiency ecc and others did not form globular clusters at all, the distribution would 
shift towards higher redshift. For example, filled circles show the mass function in the 
case when 10% of the progenitor halos form clusters with the efficiency 10 egc- 



as Mmax/ecc- For a cosmic baryon fraction, /& = Qb / &o ~ 0.13, this corre- 
sponds to the virial mass of the progenitor M v ; r = M max / (f b ecc) ~ 3000 M max . 

Figure [j] shows the distribution of the progenitor halos of the Milky Way 
calculated using the extended Press-Schechter formalism p2[ . At redshift z = 7, 
the distribution fits accurately the high end of the LF of the Galactic clusters 
renormalized by {fb^cc) 1 - In this model, the most massive cluster originates 
in the progenitor with the mass M V - 1T « 6 x 10 9 Mq and circular velocity v c w 25 
km s _1 . Correspondingly, smaller progenitors host smaller globular clusters. 

Dwarf galaxies in the Local Group provide evidence that globular clusters 
can form in such low potential wells. Even the Fornax dwarf spheroidal galaxy 
with the mass of only 2 x 10 7 M Q contains 5 clusters, while the LMC is forming 
star clusters now. Perhaps the most striking example is Li Cen, the most massive 
Galactic globular cluster. Unlike all other globular clusters, it displays multiple 
stellar populations with a significant spread of age (3 — 5 Gyr) and metallicity 
(5[Fe/H] ~ 1). Self-enrichment is unlikely, as this cluster is not special among 
others in its age and the ability to retain its own stellar winds M. If the oldest 
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Fig. 2. ART simulation of the disk galaxy formation at z = 4 (by A. Kravtsov). 
Upper panels show the gas column density and the temperature in a box of 1 h~ x Mpc 
comoving. The lower panels show the central 125 h~ x kpc around the cold disk. 

stars in to Cen had formed at its present position close to the center of the Galaxy 
(the pericenter is 1.2 kpc), the cluster would be stripped of any accumulated gas 
every 6 x 10 7 yr and would be unable to form a second generation of stars. 
The only viable scenario is that it formed within a progenitor galaxy at a large 
distance and was brought towards the center on a very eccentric orbit. 



3 Resolving globular clusters in cosmological simulations 

The complexity of gas dynamics in the high-redshift progenitors can only be 
probed with direct numerical simulations. They must involve a large range of 
scales, from the whole galaxy to the individual clusters, and include proper 
tidal interactions among the progenitors. To this end, in collaboration with A. 
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Fig. 3. ART simulation: the density, the mass accreted between successive time steps, 
and the metallicity of gas in the central 20 pc of the main progenitor. 

Kravtsov, we are analyzing the results of a very high-resolution simulation of 
the formation of a Milky Way-type galaxy. The simulation is done using the 
Adaptive Refinement Tree (ART) code Jll| , which can refine the resolution grid 
in the high-density regions by a factor 2 9 . 

Figure || shows the distribution of gas centered on the most massive progeni- 
tor at z = 4. The gas flows inward along the filaments and is shock-heated when 
it enters the dense virialized halo. A big red corona around the center in the up- 
per right panel is the hot gas at the virial temperature. The lower panels show 
a thin disk which is actively forming stars in the middle; at z = 4 its physical 
size is 2 — 4 kpc. 

In the central part of the disk the gas is not allowed to cool below 300 
K, otherwise the density becomes too high and the simulation halts. This is a 
numerical limitation of the run which can be overcome in the future. Even with 
that criterion, the density in the central resolution cell exceeds 10 5 cm~ 3 . This 
is close to the density of molecular clouds and such an object is very likely to 
form a dense star cluster. How do we determine when a globular cluster forms? 

Figure || shows the growth of density in the central cell, which has several 
noticeable peaks associated with the fast episodes of accretion. Note on the lower 
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right panel of Fig. g that the filaments stretching all the way to the center bring 
in some gas that is so dense that even if shocked, it quickly cools back. Thus 
the central region is directly accreting cold, metal-poor gas. Such episodes of 
accretion or merging are apparent in the middle panel of Fig. which shows 
the mass added between successive time steps. Some of the gas is ejected from 
the center by the hydrodynamical flows, so the mass difference can be negative. 

The two episodes at z ss 7.3 and z « 5.9 are most likely to form massive 
star clusters. The accreted gas suddenly increases the external pressure on the 
molecular clouds and triggers their collapse. The time for cluster formation is 
short, after 30 — 50 Myr the remaining gas may be thrown out by a new merger. 
After redshift of 5 the conditions become less favorable - the accretion slows 
down and the accumulated stellar mass starts to affect gas dynamics. 

The bottom panel of Fig. || shows the evolution of metallicity due to the SN 
type II ejecta. It quickly rises to 10% of solar and then varies slowly, either up 
or down. Note that the sudden events of accretion of the cold, unprocessed gas 
lowers the mean metallicity even in the very central region. If a series of globular 
clusters forms between z = 8 and 5, the younger ones are not necessarily more 
metal-rich than the older ones. Thus we should not expect a clear age-metallicity 
correlation for the clusters formed at high redshift. 

We will continue the investigation of the gas dynamics and star formation in 
all progenitor halos, with the goal of reproducing the present ages, metallicity, 
and spatial distribution of the globular clusters. 
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